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We demonstrate a novel cross-sectional deformation, 
called the radial corrugation, of multi-walled carbon nan- 
otubes (MWNTs) under hydrostatic pressure. Theoreti- 
cal analyses based on the continuum elastic approxima- 
tion have revealed that MWNTs consisting of more than 
ten concentric walls undergo elastic deformations at crit- 
ical pressure Pc — 1 GPa, above which the circular shape 
of the cross section becomes radially corrugated. Vari- 
ous corrugation modes have been observed by tuning the 
innermost tube diameter and the number of constituent 
walls, which is a direct consequence of the core-shell 
structure of MWNTs. 




Cross-sectional views of multi-walled carbon nanotube under 
high hydrostatic pressure: elliptic deformation with the mode 
index n = 2 (left), and radial corrugations with n = 5 (center) 
and n — 6 (right). The index n indicates the circumferential 
wave number of the deformed cross section. 
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1 Introduction An important mechanical feature of 
carbon nanotubes is their high flexibility in the radial di- 
rection. Radial stiffness of an isolated carbon nanotube is 
much less than axial stiffness [1], which result in an elas- 
tic deformation of the cross section on applying a hydro- 
static pressure. Such a pressure-induced radial deformation 
yields significant changes in electronic ||2]|3H!36]|7][8l|9l 
and optical 110, , 11 , 12. 13| properties, indicating the rele- 
vance of the deformation in carbon nanotube applications. 
Thus far, many experimental and theoretical studies have 
been carried out on this issue lfT4 15 16 17 18 19 2(1211 
|22l|23]|24l|25]|26, 27, 28, 29, 30, 31,32,33,341. Most of them 
focused on single-walled nanotubes (S WNTs) and their bun- 
dles, and revealed flattening and polygonalization in the 
cross section of S WNTs under pressures of the order a few 
GPa 1101161. Contrary to the intensive studies on SWNTs, 
radial deformation of multiwalled nanotubes (MWNTs) still 
remains to be explored. Intuitively, the multiple-shell struc- 
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ture of MWNTs is thought to enhance the radial stiffness of 
MWNTs. However, when the number of concentric walls is 
much greater than unity, outside walls have large diameters 
so that external pressure may lead to a mechanical instabil- 
ity in the outer walls. This local instability implies a novel 
cross-sectional shape transition of MWNTs different from 
the cases of SWNTs. 

This paper shows a new class of radial deformation, 
called radial corrugation, of MWNTs consisting of more 
than ten concentric walls. In a corrugation mode, outside 
walls exhibit wavy structures along the circumference as 
depicted in Titlefigure. We demonstrate that various cor- 
rugation modes can take place above critical hydrostatic 
pressure Pc{N, D), in which pc depends on the innermost 
tube diameter D and the number of constituent walls N . 
It should be emphasized that radial corrugation we have 
found is a direct consequence of the core-shell structure 
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Figu re 1 Illustrations of geometric parameters of the continuum elastic-shell model, (a) Sketch of cross section of MWNT 
consisting of N cylindrical walls with thickness h. Each wall is labeled by the index i. (b) Enlarged view of a partial cross 
section of a cylindrical wall before deformation. The circumferential line element and the centroidal circle C are 
depicted, (c) Enlarged view of a partial cross section of a deformed wall under high pressure. The line element is elongated 
in the circumferential direction, and the normal to the centroidal curve rotates with an angle /3. 



of MWNTs, and thus is differernt inherently from simple 
radial collapse observed in SWNTs. 

2 Continuum elastic-shell theory 

2.1 Outline The stable cross-sectional shape of a MWNT 
under a hydrostatic pressure p is evaluated by using the 
continuum elastic theory for cylindrical shells i35': 36l[37l 
[38l . The mechanical energy [/ of a MWNT per unit axial 
length is written as 



u = u\v, u,{p, e), v,{p, 0)] = UD + Ui + n, 



(1) 



where Ud is the deformation energy of all concentric walls, 
Ui is the interaction energy of all adjacent pairs of walls, 
and n is the potential energy of the applied pressure. All 
the three energy terms are functions of p and the defor- 
mation amplitudes Ui{p,6) and Vi{p,6) that describe the 
radial and circumferential displacements, respectively, of 
the ith wall. See Eq. ^ below for the precise definitions 
of Ui and Vi . 

Our objective is the optimal displacements Ui and Vi 
that minimize the total energy U under a given pressure p. 
To this aim, we apply the variation method to U with re- 
spect to Ui and Vi, and then obtain the stable cross-sectional 
shape of a MWNT under p ||39l . This strategy requires to 
derive explicit forms of Ud, Uj and f2 as functions of Ui, 
Vi and p, which we shall resolve in the subsequent discus- 
sions. 

2.2 Strain-displacement relation We first consider 
the circumferential strain of a hollow cylindrical shell due 
to cross-sectional deformation. Suppose a circumferential 



line element of length di lying within the cross section of 
the shell with thickness h (See Fig. [T]i. The tilde ( ~ ) at- 
tached to £ means that we consider the quantity at arbitrary 
point within the cross section. We shall see later {i.e., in 
Eq. (O ) that the strain at arbitrary point is determined ap- 
proximately by the strain just on the controidal circle de- 
noted by C in Fig.[T| This fact allows us to yield a simplified 
relation between the strain and displacements of the shell 
as given m Eqs. dill and 

The extentional strain e of the circumferential line ele- 
ment is defined by 



d£* - d£ 

d! ■ 



(2) 



Here d£ = fdO, and d£* is the length of the line element af- 
ter deformation (The asterisk symbolizes the quantity after 
deformation). The coordinates x* , y* of the element after 
deformation is given by 



x*{e) = [f + u{e)]cose ~d{d)sm9, 
y* (9) = \r + u{e)] sin 9 + i{9) cos 9, 



(3) 



where u and v are components of the displacement vector 
in the radial and circumferential directions, respectively. It 
thus follows that 



d£*] ^{dx*f + {dy*f 



f"^ + 2f [u + v) + {u + vf + [u - vf {d9) 



(4) 
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FromEq. (|2]i, we have e+1 ~ d£* /d£. Hence, squaring 
the both sides and then rearranging the resuh give 



2 2 



d£ 



(5) 



For e <C 1, the term e"^ can be omitted. Hence, we have 
from Eqs. (01 and ^ that 











f u + v' 















(6) 



where u' = du/dd, etc. The last term in Eq. @ is as- 
sociated with the rotation of the Une element due to de- 
formation. The rotation angle j3 consists of two parts: i) a 
clockwise component du/d£ = du/{fd9) due to the spa- 
tial variation of u in the circumferential direction, and ii) 
a counterclockwise one v/f due to the circumferential dis- 
placement of the element. Combination of the two parts 
gives 



(7) 



which has a positive value in the counterclockwise sense. 

The formula (|6]l is valid for arbitrary large rotation (3. 
Particularly when e and (3 are both sufficiently small (but 
finite), we may neglect the second term in the right side in 
Eq. (|6]l. (Here we exclude the possibility that |u| or \v'\ is of 
the order of f or larger.) We further assume that normals to 
the undeformed centroidal circle C remain straight, normal, 
and inextensional during the deformation (See Fig. [T]). As 
a result, u and v are expressed by 



u — u and v — v + z(3. 



(8) 



where u and v denote the displacements of a point just on 
C, and z is a radial coordinate measured from C. By sub- 
stituting (|8]l into (|6]l, we attain the strain-displacement re- 
lation such as l l40l 

(9) 



e{z,e) = e{0) + zK{9), 
with the definitions: 

u + v' 1 I u' — V 



and K 



(10) 



Here r is the radius of the undeformed circle C. The results 
(|9]l and ( fTOl i state that the circumferential strain at arbitrary 
point in the cross section is determined by the displace- 
ments u{6) and v{9) of a point just on the undeformed 
controidal circle C. 

2.3 Deformation energy We are ready to derive the 
explicit form of the deformation energy Ud- Suppose the 
ith cylindrical wall of a long and thin circular tube with 
thickness h. A surface element of the cross-sectional area 
of the wall is expressed by ridOdz. The stiffness k of the 
surface element for stretching along the circumferential di- 
rection is given by 

fc = — (11) 

1 — 



where E and v are Young's modulus and Poisson's ratio, 

respectively, of the wall. Thus, the deformation energy Uj^ 
of the ith wall per unit axial length is written as 

, fh/2 .27r 

/ i,{z,efdzde. (12) 

h/2 Jo 

From (fT2l i and (|9]l, we obtain 

khr, f^'" , ,„ kh^r '"^'^ 



" 2 



U 



efde 



24 



Kfde, 



(13) 



which tells us the dependence of Ud ~ J2iLi U^'^ on 
u,{e) andv,{e). 

2.4 Inter-wall coupling energy Secondly, the ex- 
plicit form of Ui is described by 







N 

E- 

1=2 



27r 



{Ui - dO. 



(14) 



The vdW interaction coefficients c,, are functions of r, and 



Tj and defined by 11381 

/100l7r£cri2 



-F, 



13 



11207r£crS 



(15) 



where a denotes the chemical bond length between neign- 
bouring carbon atoms within a layer (a = 0.142 nm), and 
£ and a are the parameters that determine the vdW interac- 
tion between two layers {e = 2.968 meV and a = 0.3407 
nm) j44 1. In Eq. ( fTSl l. we have set 

r.7r/2 







[ 1 - K,j cos^ 61) 



(16) 



ij — Arirj/{ri + rj) , and the wall spacing |r,; 



with K^j = 4nrj/{ri + r^f 

ri±i \ = 0.344 + 0.1c--°/2 nm according to Ref. 1421 . 

2.5 Pressure-induced energy We finally derive an 
explicit form of i7. Since i7 is the negative of the work 
done by the external pressure during cross-sectional defor- 
mation, it is expressed as 

fi = -pinrj, ~ S*). (17) 
Here, S* is the area surrounded by the A'^th wall after de- 
formation (the sign of p is assumed to be positive inward). 
It is obtained by the Une integral 5* = ^ x dr^\, 

or equivalently 



S* = 



1 



2tt 



{xWn - V*nX*n') de. (18) 

From Eqs. (fTsT l and (O as well as the periodicity relation 

J^^ v^'dO = 0, we obtain 



rNUN 



uIj + v% - un'vn + unvn' 



de. 



Table 1 Values of elasticity parameters used in this article. 

Young's modulus Poisson's ratio Wall thickness 



E = 1 [TPa] 



u = 0.27 



h = 0.34 [nm] 
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Number of walls: N Innermost tube diameter: D [nm] 



Figure 2 Left: Wall-number dependence of the critical pressure Pc{N) for various innermost tube diameters D. Each 
kink appearing at = A^o = 17 for D = 4 nm and Nq = 30 for D = 5 nm separates the elliptic deformation phase 
{N < Nq) and the radial corrugation phase {N > Nq). For D > 6 nm, only the elliptic deformation occurs independently 
of N. Right: The D-dependence of Pc{D) for various Ns. Radial corrugation occurs at £> < Dq, where Z?o indicates the 
position of the kink appearing for N > 20. 



3 Critical pressure curve 

3.1 Evaluating critical pressure pc Our aim is to 

evaluate the critical pressure pc above which the circu- 
lar cross section of a MWNT is elastically deformed into 
non-circular one. This is achieved by the following pro- 
cedure 1391 . First, we decompose the radial displacement 

terms as Ui(j>,9) = uf^'^ip) + 5ui{6), where uf'^p) in- 
dicates a uniform radial contraction at p < pc and 5ui{6) 
describes a deformed, non-circular cross section observed 
just above Similarly, we can write Vi{p,9) = 6vi{9), 
since vf'\p) = at p < pc- Next, we apply the variation 
method to U with respect to Ui and Vi, which results in a 
system of 2N linear differential equations with respect to 
6ui and 5vi. Then, substituting Fourier series expansions 
of 5ui and 5vi into the differential equations, we obtain the 
matrix equation Cu = Q. Here, the vector u consists of 
the Fourier components of 5ui and 5vi with all possible 
i and all Fourier components labeled by n, and the ma- 
trix C involves one variable p as well as elasticity param- 
eters. Finally, solving the secular equation det(C) = 
with respect to p, we obtain a sequence of discrete values 
of p; among these ps, the minimum one serves as the criti- 
cal pressure pc immediately above which the circular cross 
section of MWNTs becomes radially deformed. 

3.2 Wall-number dependence Figure|2]shows nu- 
merical results of critical pressures pc- In the left-side fig- 
ure, Pc is plotted as a function of the number of concentric 
walls N with the innermost tube diameter D being fixed. 



For all D, pc increase with N followed by a slow decay or 
saturating constant value. The increase in pc in the region 
of small is attributed to the enhancement of radial stiff- 
ness of the entire MWNT by encapsulation. In contrast, the 
decay in pc indicates local instability in outside walls that 
causes radial corrugation of MWNTs with ^ 1. 

To be noteworthy is the occurrence of a kink mpcXN) 
for small D. In the curve of Pc{N) for D = 4 nm, for 
instance, a kink appears at A^ = A'o = 17 to the right of 
which the value of pc decay monotonically. The most strik- 
ing is the fact that this kink determines the phase boundary 
between the elliptic deformation phase and the radial cor- 
rugation phases (as indicated by red arrows in the left-sided 
figure of Fig. 2). We have found that A^-walled nanotubes 
with A^ > A'o exhibit some radial corrugation just above 
Pc, while those with A^ < A'o show only the elliptic de- 
formation. This senario also holds for the case of Z? = 5 
nm, where the phase boundary kink occurs at A'o = 30. 
Which kind of corrugation mode (i.e., what number of the 
mode index n) takes place is dependent on the values of 
N{> No) and D, as clarified later in Fig. [3] 

3.3 Innermost-tube-diameter dependence The 
right-side figure in Fig. |2] illustrates the D-dependence of 
Pc{D) for various Ns. Similarly to the left-side figure, a 
kink appears for each curve except for A^ = 10. This kink 
again serves as the phase boundary between the elliptic and 
corrugation phases; if a kink occurs at Dq, then it separates 
radial corrugation phases {D < Dq) from the elliptic de- 
formation phase (D > Dq). 



Copyright line will be provided by the publisher 



pss header will be provided by the publisher 



5 




Number of walls: N 



Figure 3 Left: Three-dimensional plot of the critical pressure Pc{N, D) in the N-D plane. A ridge line extending from 
the top to the bottom of the surface corresponds to the phase-boundary kinks that separate the elliptic phase from multiple 
corrugation phases. Right: Phase diagram of the pressure-induced cross-sectional deformation of MWNTs. 



It is interesting that in the elliptic phase, all the results 
(except for TV = 10) collapse onto a single curve. This 
means that the number of walls N is irrelevant in determin- 
ing pc for MWNTs with large D. In other words, encapsu- 
lation by increasing N gives no contribution to the radial 
stiffness of those MWNTs with large D. This fact implies 
the availiability of simplified equations for estimating pc of 
MWNTs with large D; establishing such simplified equa- 
tions is currently in progress. 

4 Radial corrugation of lUIWNTs All the numeri- 
cal results of pc are summarized into the three-dimensional 
plot Pc{N, D); see the left-side in Fig. [3] This plot makes 
clear the optimal values of N and D that maximize (or 
minimize) the radial stiffness of MWNTs. We also obverve 
a ridge line extending from the top to the bottom of the 
surface. The ridge line corresponds to the phase-boundary 
kinks discussed in the previous section. In fact, multiple 
corrugation modes take place in the region to the left of 
the ridge line, as demonstrated in the phase diagram; see 
the right-side of Fig. [3] Interestingly, multiple corrugation 
modes are formed depending on the values of N and D. We 
see that smaller D and larger N favor corrugation modes, 
in which larger N yield higher corrugation modes with 
larger n. 

5 Concluding remarks We have demonstrated the 
presence of multiple radial corrugations peculiar to MWNTs 



under hydrostatic pressures. Theoretical investigations based 
on the continuum elastic theory have revealed that MWNTs 
consisting of a large number of concentric walls undergo 
elastic deformations at critical pressure Pc ^ ^ GPa, above 
which the cross-sectional circular shape becomes radially 
corrugated. A phase diagram has been established to ob- 
tain the requisite values of N and D for observing a de- 
sired corrugation mode. It is hoped that our results should 
be verified by high-pressure experiments on MWNTs as 
well as atomic-scale large-scale simulations II43II44||45]| . 

Another interesting subject is to explore the structural 
deformation effect on the quantum-mechanical properties 
of TT electrons moving in the corrugated carbon walls. It 
has been known II46II47B8I that mobile electrons whose 
motion is confined to a two-dimensional curved thin layer 
behave differently from those on a conventional flat plane; 
this results from the occurrence of effective scalar- 1491 
l50l and vector- lISTI potential energies induced by geomet- 
ric curvature of the underlying layer. Hence, quantum na- 
ture of electrons moving in the corrugated nanotube will 
be strongly affected by geometric curvature of the walls, 
which remains thus far unsettled. Intensive studies on the 
issues mentioned above will shed light on novel MWNT 
applications based on cross-sectional deformation. 
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